Particle systems with quasi-homogeneous initial states and 
their occupation time fluctuations* 

Tomasz Bojdecki 1 Luis G. Gorostiza^ Anna Talarczyk 1 

tobojd@mimuw.edu.pl lgorosti@math.cinvestav.mx annatal@mimuw.edu.pl 

1 

Abstract 

Occupation time fluctuation limits of particle systems in WL d with independent motions (sym- 
metric stable Levy process, with or without critical branching) have been studied assuming initial 
distributions given by Poisson random measures (homogeneous and some inhomogeneous cases). In 
this paper, with d = 1 for simplicity, we extend previous results to a wide class of initial measures 
obeying a quasi-homogeneity property, which includes as special cases homogeneous Poisson mea- 
sures and many deterministic measures (simple example: one atom at each point of Z), by means 
of a new unified approach. In previous papers, in the homogeneous Poisson case, for the branch- 
ing system in "low" dimensions, the limit was characterized by a long-range dependent Gaussian 
process called sub- fractional Brownian motion (sub-fBm), and this effect was attributed to the 
branching because it had appeared only in that case. An unexpected finding in this paper is that 
sub-fBm is more prevalent than previously thought. Namely, it is a natural ingredient of the limit 
process in the non-branching case (for "low" dimension), as well. On the other hand, fractional 
Brownian motion is not only related to systems in equilibrium (e.g., non-branching system with 
initial homogeneous Poisson measure) , but it also appears here for a wider class of initial measures 
of quasi-homogeneous type. 
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1 Introduction 



In a series of papers [3l[31 [5j [61 [Tl [HI [10] we studied particle systems in R d starting from a configuration 
determined by a random point measure u, and independently moving according to a standard a-stable 
Levy process (0 < a < 2). In some models the particles additionally undergo critical branching. The 
evolution of the system is described by the empirical process N = (Nt)t>o, where Nt(A) is the number 
of particles in the set A C M. d at time t. The main object of interest is the limit of the time-rescaled 
and normalized occupation time fluctuation process Xt defined by 

1 f Tt 

X T (t) = — (N s - EN s )ds, t > 0, (1.1) 

b T Jo 

as T — > oo (i.e., as time is accelerated), where Ft is a suitable deterministic norming. The process Xt is 
signed measure- valued, but we regard it as a process with values in the space of tempered distributions 
S'(M. d ) for technical convenience, and also because in some cases the limit is genuinely 5'(IR d )-valued. 
In all the cases considered in the abovementioned papers the initial measure v was a Poisson field, 
homogeneous or not. This assumption permitted to investigate convergences conveniently with the 
help of the Laplace transform (due to infinite divisibility) . The results always exhibited the same type 
of phase transition: for "low" dimensions d the limit process was the Lebesgue measure multiplied by 
a real long-range dependent process, whereas for "high" dimensions the limit was an 5'(M d )-valued 
process with independent increments. 

A natural question is what happens for non-Poisson initial measures v. Milos 127J[28|[29] considered 
(critical) branching systems where v was an equilibrium measure (see [IS])- In that model the limits 
have a similar dimension phase transition; moreover, for high dimensions they are the same as in the 
homogeneous Poisson case, while for low dimensions they are different. The conclusion was that in 
low dimensions the occupation time fluctuation process "remembers" the initial state of the system. 
Since the equilibrium states of the branching system are somewhat similar to homogeneous Poisson 
measures (they are infinitely divisible random point measures with uniform intensity; distributions of 
this kind are called "equilibrium distributions of Poisson type" in [26]), the Laplace transform method 
was also useful in [27], [28| [29] . 

The aim of the present paper is to investigate what happens with initial measures of other types, 
for example, some measures that are deterministic or almost deterministic. For simplicity we consider 
d = 1 and assume that the motions are either without branching or with the simplest critical binary 
branching. In [2J [5] we proved for such motions, with general d, that if v is a homogeneous Poisson 
measure, then the following results hold (where A denotes Lebesgue measure and K is a different 
constant in each case): 

in the non-branching system: 

if d < a, then Xt converges in distribution (in C([0, r], 5'(IR rf )) for any r > 0) to a process K\£, 
where £ is a fractional Brownian motion; 

if d = a, then the limit process is KX/3, where /3 is a standard Brownian motion; 
if d > a, then the limit is a time-homogeneous Wiener process in <S'(R ); 
in the branching system: 

if a < d < 2a, then the limit is KX^, where £ is a sub- fractional Brownian motion (the case d < a 
requires a slightly different treatment based on high-density models, see [9]); 
if d = 2a, then the limit is -ff A/3; 

if d > 2a, then the limit is a time-homogeneous Wiener process in 5'(IR d ), different from the one 
in the non-branching case. 
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In this paper we define a class Ai of initial measures v which contains in particular homoge- 
neous Poisson measures (which are "completely random" [23]), and quasi- homogeneous deterministic 
measures (e.g., the measure defined by one atom at each j G Z, which is "completely determinis- 
tic"), and we develop a unified approach that permits to obtain limits of Xt for all v G Ai. By a 
quasi-homogeneous deterministic point measure on R we mean any measure defined by the following 
procedure: Given a positive integer k, in each interval [j,j + 1), j G Z, we fix k points. For a general 
v G Ai, each interval [j,j + 1) contains 6j points chosen at random, and 9j, j G Z, are i.i.d. random 
variables (see Section 2 for a rigorous definition). The main feature of those measures is this form of 
quasi-homogeneity and independence on the family of intervals [j,j + 1). 

For each v G Ai we obtain the limit of the corresponding Xt and in this way we recover the 
results of IU [5] for the homogeneous Poisson case (for d = 1 , but there is no doubt that the results 
for higher dimensions can be obtained analogously), and we also derive limits for many other initial 
measures. It seems interesting that the idea of the proofs in this general framework is simpler than 
that in our previous papers, and is based on the central limit theorem. This is a significant change of 
methodology. However, some technical points in those papers are employed again here. An argument 
using the non-linear equation associated with the occupation time of the branching system, which 
can be obtained by means of the Feynman-Kac theorem, again plays an important role, but now in 
a different way: it is a key step in moment estimates in order to apply the Lyapunov theorem in the 
branching case. The equilibrium measures for the branching system do not belong to Ai because the 
branching introduces spatial dependence. 

Some of the results we obtain are unexpected. It turns out that the only case where new limits 
appear is the non-branching case with (d =)1 < a. They have the form KA£, where £ is the sum of two 
independent processes, one of them is a sub-fractional Brownian motion (see (12. 3p ). and the second one 
is a new (centered continuous with long range dependence) Gaussian process (see (|2.4p ). The process 
£ depends on the initial measure v only through E9q and Var#o- In particular, for a deterministic 
initial measure this process reduces to a sub-fractional Brownian motion, and in the homogeneous 
Poisson case (as well as for any u with E9q = Var#o) it yields a fractional Brownian motion (see 
Theorem 12. 2p . This result seems surprising since in all earlier papers sub-fractional Brownian motion 
was related only to branching systems, and was consequently attributed to the branching, but now, 
in the present context, this process turns out to be more "natural" than fractional Brownian motion. 
On the other hand, fractional Brownian motion, which is typically related to systems in equilibrium 
(in particular the non-branching system with initial homogeneous Poisson measure), now appears also 
for a wider class of quasi-homogeneous initial measures, as noted above. 

In all the remaining cases the limits are (up to constants) the same, and with the same normings 
Ft, as those recalled above for homogeneous Poisson models. 

The results show that within the class Ai the fluctuations caused by the branching are so large 
that Xt "forgets" the randomness of the initial state of the system (it "remembers" EOq only). On 
the other hand, for low dimensions it does distinguish between v G Ai and the equilibrium initial state 
(which is not in Ai). Another conclusion is that for high dimensions (which for d = 1 amounts to small 
a), the fluctuation process "forgets" the initial measure, as long as it is in some sense homogeneous 
(i.e., i/fM), and this property holds for branching and non-branching systems; it is also preserved 
for branching systems in equilibrium. 

In this paper we are interested mainly in identifying the limit processes, therefore we have not 
attempted to prove convergences in their strongest, functional form; in most cases we prove only 
convergence of finite-dimensional distributions. Presumably, convergence in distribution also holds in 
C([0, r], <S'(R)) for any r > 0. As an example, we give one result of this type (Proposition 12. 6p . 

Some other papers related to occupation times of particle systems and superprocesses are [21 HH 
H2| [T3"l dH [151 El 1211 1221 1231 [Ml EH] . For example, occupation time limits for branching random walks 
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on the d-dimensional lattice are discussed in [2]. 

We have not found results in the literature concerning occupation times for particle systems starting 
from a deterministic or quasi-deterministic point measure. Some kinds of quasi-homogeneity of initial 
configurations for systems of independent particles, different from those in this paper, appear in other 
contexts in [33j and [20j (see Remark 12.51 (f)). It may be that systems of independent particles with 
a-stable motion and the initial conditions of [33J lead to the same results as with initial homogeneous 
Poisson distribution. 

The following notation is used in the paper. 

<S(R): space of C°° rapidly decreasing function on R. 

S'(R): space of tempered distributions (topological dual of <S(R)). 

( , ): duality on S'{R) xS(R). 

=>f. weak convergence of finite-dimensional distributions of 5'(R)-valued processes. 
Pt(x): transition probability density of the standard symmetric a-stable Levy process. 
%: semigroup determined by pt, i.e., Ttf =pt*<p- 
G: potential operator determined by pt for a < 1, i.e., 

G<p{x) = r°rMx)db = C a I , ^ a dy, (1.2) 
Jo Jr \x - y\ L a 

where 

a 2 a 7T 1 /2r(f)' 

Generic constants are written C, Cj, with possible dependencies in parenthesis. 

In Section 2 we describe the particle system, formulate the results and discuss them 
contains the proofs. 



2 Results 

We start with detailed description of the particle system. 

Let 9 be a non-negative integer-valued random variable with distribution 

P{6 = k) =p k , k = 0,1,2,..., (2.1) 

such that E8 3 < oo. This moment condition is a technical assumption satisfied by all cases of interest 
in this paper, but we suppose that finiteness of the second moment could be sufficient. 

Let 9j,j G Z, be independent copies of 9, and for each j G Z and k = 1,2,..., let f? k = 
(p J k v . . . , k ) be a random vector with values in [j,j + l) fc . We assume that (9j, (pl)k=i,2, ■ ■ ■)> 3 ^ Z, 
are independent. These objects determine a random point measure v on R in the following way: For 
each j, 9j is the number of points in the interval [j,j + 1), and for each k, if 9j = k, the positions of 
those points are determined by pi. In other words, 

0j 

^ = EE^,n> (2-2) 

jgZn=l 

where Kj n = p@ and 5 a is the Dirac measure at a G R. We denote by Ai the class of all such 
measures v. 
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Remark 2.1 (a) If 9 = k and, for each j, the f? k are not random, then v is a quasi-homogeneous 
deterministic measure mentioned in Introduction. The simplest example is v = Yljez^j- 

(b) If 6 is a standard Poisson random variable and, for each j, fP k -, , . . . , fP k k are independent, 
uniformly distributed on [j,j + 1), then v given by l\2.2\) is the homogeneous Poisson point measure 
(with intensity measure A). 

Fix a £ (0, 2] and v E A4. Assume that at the initial time t = there is a collection of particles 
in M. with positions determined by v. As time evolves, these particles move independently according 
to the standard a-stable Levy process. We consider systems either without branching, or with critical 
binary branching (i.e., or 2 particles with probability 1/2 each case) at rate V. For the corresponding 
empirical process, we define an 5'(R)-valued process Xt by (jl.ip . 

Before stating the first theorem we recall the definition of sub-fractional Brownian motion. A 
sub-fractional Brownian motion with parameter H (0 < H < 1) is a centered continuous Gaussian 
process ( H with covariance 

C H (s, t) = E(?( t H = s 2H + t 2H ~\{{s + t) 2H + \s- t\ 2H ), s,t>0. (2.3) 

See [21 [35] for properties of this process. It appears in [16] in a different context, and it has also been 
investigated in [H ED [361 EH EE] . 

We will need another centered Gaussian process $ H with covariance 

Q H (s, t) = isgn(2tf - 1)(( S + t) 2H - s 2H - t 2H ), s,t>0, (2.4) 

(0 < H < 1). Existence of this process for H ^ 1/2 follows from the formula 

Q H (s,t) = C f [\r + r') 2H - 2 dr'dr = C 1 f f \ e^ 1 ^ e^ 1 ^ dxdr' dr, 
Jo Jo Jo Jo Jr 

which implies positive-definiteness of Q H . 

Theorem 2.2 For the system without branching, 
(a) if 1 < a and 

F T = T 1 - 1 / 2 ", (2.5) 

then 

X T ^f Ki\{\fEeC, H + War 9$ H ), (2.6) 
where £ , i? are independent, H = 1 — l/2a, and 



(b) if a = 1 and 



( r(2 - 2H) 
\2iraH(2H - 1) 



1/2 



F T = v/TlogT, (2.7) 
then 

where (3 is a standard Brownian motion in M., and 



A 2 - 
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(c) if 1 > a and 

F T = VT, (2i 

then 

X T ^ f X, (2.£ 
where X is an S' (R) -valued homogeneous Wiener process with covariance 



E{X(t),ip){X(s),i/>) = 2E9{s At) / ip(x)G^(x)dx, <p,il> € S(NL). (2.10) 

Jr 

Remark 2.3 (a) Let 

£ H = VE9( H + \/Var0tf H (2.11) 

be the process in Theorem 2.2(a). If E9 = Var#, in particular if v is homogeneous Poisson (see Remark 
2.1(b)), then is, up to a constant, a fractional Brownian motion with Hurst parameter H, i.e., it 
has covariance C{s 2H + t 2H -\s-t\ 2H ). Thus we recover Theorem 2.1 of [3]. On the other hand, if 
9 is deterministic, then £^ is a sub-fractional Brownian motion. Moreover, in general randomness of 
the p's in the definition of v € Ai (see (2.2)) does not play any role in the limit. 

(b) The long time dependent behavior of Gaussian processes is usually characterized by the co- 
variance of increments of the process on intervals separated by distance r, as r — > oo. For the process 
•d H that behavior is asymptotic decay like t 2H ~ 2 (the same as for fractional Brownian motion), for 
sub-fBm it is r 2H ~ 3 (see [3]). So, the long time dependent behavior of the process £ H in Theorem 
12.2( a) is determined by d H in all cases where 9 is random. 



Theorem 2.4 For the branching system, 
(a) if 1/2 < a < 1 and 

p T = r (3-l/a)/2 5 ( 2 _ 12 ) 



Xt =>f K 3 X( H , 



then 

where C, H is a sub-fractional Brownian motion with parameter H = (3 — 1 /a) jl, and 

K 3 



/E9VT(2 - 2H)\ 1/2 
\2iraH(2H - 1) J 



(b) if a = 1/2 and Ft = \/T logT, then 

X T =>/ -fC 4 A/3, 
where (3 is a standard Brownian motion, and 



2V 

K 4 = ^j-EO; 

(c) if a < 1/2 and Ft = VT, then 

Xt =^f X, 

where X is an S' (R) -valued homogeneous Wiener process with covariance 

E{X(t),<p){X(s),i/>) = E9(s At) / (2tp(x)Gil>(x) + V{G<p(x)){G^{x)))dx, <p,ip G <S(R). (2.13) 
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Remark 2.5 (a) In the branching case the results are, up to the constant E6 in the limits, the same 
as in the homogeneous Poisson case (Theorems 2.2 in [3] and [5]). 

(b) The condition a < 1 in part (a) of the last theorem corresponds to a < d in [3] and |5]. In 
the homogeneous Poisson case for d < a, we obtained limits of the same form as for a < d < 2a, by 
introducing high density, i.e., considering systems with initial intensity Ht\, Hj- — > oo sufficiently fast 
[9]; the high density counteracts the tendency to local extinction caused by the critical branching. The 
same procedure can be applied in the present case, yielding the limits for 1 < a < 2, if the intervals 

+ 1) are replaced by [j/H T , (j + l)/H T ). 

(c) As in [4 J and [5j, Theorems 2.2 and 2.4 can be extended to systems in M , where the intervals 
[j,j + 1) are replaced by cubes [j, j + l) d . 

(d) Comparing parts (a) of Theorems 2.2 and 2.4, we see that the branching weakens the influence 
of the initial configuration. 

(e) The previous results show that sub-fractional Brownian motion is a "natural" process for our 
model. So far it had appeared only in the context of branching systems, but now we see that it 
is intrinsically related to the non-branching systems as well for a large class of initial conditions. 
Fractional Brownian motion occured before only in the case of systems in equilibrium, but now it also 
appears wherever E6 = Var#. 

(f) Theorems 2.2 and 2.4 can also be extended to other models. For example, in [20] a model is 
studied in a different context with independent a-stable motions without branching and initial posi- 
tions of particles (j + p)j^z, where p is a random variable uniformly distributed on [0, 1], independent 
of the motions. It is easy to see, by a standard conditioning argument (considering the characteris- 
tic function and conditioning on p), that for models with or without branching and with this initial 
configuration, the limits are the same as for deterministic p, i.e., they are given by Theorems 2.2 and 
2.4. 

We have formulated Theorems 2.2 and 2.4 with convergence of finite-dimensional distributions only, 
as we are mostly interested in the limit processes, but we have no doubt that functional convergence 
also holds. As an example, let us consider the cases of large a. For simplicity we assume that the 
initial configuration is such that 8 = 1 and p\ 1 is uniformly distributed on [j, j + 6 Z. 

Proposition 2.6 For the model described above, the processes Xt in Theorems 2.2(a) and 2.4(a) 
converge in law in C([0,t],S'(M.)) for any r > 0. 

3 Proofs 

3.1 Auxiliary facts related to the stable density 

We will often use the self-similarity property of the transition density pt of the standard a-stable 
process in M: 

Pat ( x ) = a - l l a p t (a- 1/a x), xeR, a > 0. (3.1) 

Recall that 

Pl(x) < 1 + | ^ |1+Q - (3.2) 
Since pt{ ) is decreasing on IR + and symmetric, then by (|3.1j) we have 

Pt(x + y)< g t (x) := j » . f ~ 2 , * e % \v\ < 1- (3-3) 
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Denote 

<t>m(x) = 1 + | g | m . m>0. (3.4) 

For ip G <S(R) we have |y?(x)| < C((p,m)(f> m (x). This, and an obvious inequality, 

1 1 + \v\ m 

< C(m)— m > 0, (3.5) 



1 + |x + y| m 1 + \x\ 

imply 

77|^| (rc + y) < C m (a)Tt4> m {x), |y|<a,m>0, (3.6) 
In the sequel we will use various versions of this estimate, e.g., 

Tt(\<p\T s \i(>\)(x + y) < C m {a)Tt{<t>mTs4> m ){x), if,iP 6 5(1). 

We will also need the following estimate (|22j. Lemma 5.3) for the potential operator G (see 
(ll.2j) . (jl.3p ). If a < 1, q > 1, and is a measurable function on M such that 

\<p(x)\ < Ci— ^pj-, xeR, 

then 

|Gy(g)| <C 2 - , , 1 -ll _ a , xeK. (3.7) 



1 + \x\ 



3.2 Scheme of proofs 



The proofs of Theorems 2.2 and 2.4 are based on the central limit theorem and follow the scheme 
described presently. 

Let N x denote the empirical process of the system (with or without branching) started from a 
single particle at x, and N^\j E Z, be the empirical process for the particles which at time t = 
belong to [j,j + 1), i.e., 

9j 

Nii) = Y^N K ^, (3.8) 

n=l 

according to the description at the beginning of Section 2 (see (2.2)). Note that N^\ j 6 Z, are 
independent. 

The process X? defined in (jl.ip can be written as 

1 r Tt 

X ^) = Y.T ^-EN^)ds. (3.9) 

The first step in our argument is to prove that for any ip,ip G <S(M), and s,t > 0, 

Jim E{X T (t),<p){X T (8),4) = E{X(t),<p){X(8),1>), (3.10) 

where X is the corresponding limit process. Without loss of generality we may assume that <p, ip > 0. 
Using (3.9) we have 

E(X T (t),ip)(X T (s),iP) 

r-Tt f-Ts 



1 pit. plS 

Ep2/ / E{N^\ip){N^^)dr'dr 

i pTt pTs 

V-2 / / E{N(?\v)E{Ny^)dr'dr. (3.11) 

77^ T JO JO 
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Using (|3.8p . (|2.1|) and the fact that E(Nf, (p) = Ttf(x) in both non-branching and (critical) branching 
cases, and defining, for x E M, n < k, random variables 



h k ,n( x ) = Pkli ~ x > 
where [x] is the largest integer < x, we rewrite (|3.11[) as 



(3.12) 



k=0 n=l k=0 n,m=l 

OO k OO I 



(3.13) 



k=0 n=l £=0 m=l 



where 



I{T;k,n) = 



Tt r Ts 



1 



"t j o jo 



E(Nr +hk - n{x \ ip)(N* +hk -" {x \ iP)dxdr'dr, 



Tt r Ts 



"t j o jo 



E(T;k,n,m) = — tt 



"t j o Jo 



Tt r Ts 



(3.14) 



1 



Tt rTs 



E (%ip(x + hk n ( x ))'Tr''il)(x + hk, m (x))) dxdr' 'dr (3.15) 



(in the first equality for E we used independence of systems starting from different points), 



M(T;k,n;£,m) 



Tt rTs 



1 

E T jo Jo 

rTt rTs 



E ETM^jETr^ Jdr'dr 



= -=2 [ [ f E%(f(x + hk,n(x))ET r 'ip(x + hi^ m (x))dxdr' dr. (3.16) 

-f*T JO Jo Jr 

Note that 

\h k ,n{x)\ < 1, x € R. (3.17) 

In each case we will show convergence of I,E and M, thus proving (|3.10p . (It will be shown that 
/, E, EI are bounded, so the passage to the limit in each sum in (|3.13p is justified). 
Next, we show that 

(X(t),<p) => (X(t),<p), <peS(R),t>0. 
To this end, by (|3.9|) and (|3.1Up it suffices to prove that the Lyapunov condition 



1 r Tt 

lim^-^E / ((N^, ( p)-E(NjP\ ( p))dr 

icit T JO 



is satisfied, and this property will follow if we show that 

,-Tt 3 

V— / ( / 







i / rTt \ 3 
lim > -^tE / (M J \(p)dr = 0, t > 0, <p G S(W), p>0. 
T ^°°~^ F T \Jo J 



(3.18) 
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It is clear that convergence in law of linear combinations J2T=i a k{XT(tk), ¥>fc) can be obtained 
analogously from f|3. lOj) and (|3.18p . thus establishing the claimed convergence Xt j X . 

In order to give (I3,18|) a more tractable form we use (|2.ip . (13. 8p . and the trivial inequality (ai + 
. . . + a k ) 3 < 3k 2 (af + . . . + a|), a\ . . . ,a k > 0, obtaining 

r-Tt \ 3 



oo / k fTt j \ 

= EwE^ E/ <^"^>*- 

— T k=0 \n=l Jo 



< 3E8 3 sup — E[ (Nr +hk ' n{x \ip)dr) dx 



T 

n<-k 

(see (|3.12p ). So, to prove f|3. 18j) it suffices to show that 

fTt . \ 3 

SUP 4r / E( I 

T-too 



lim sup 2- [ e( [ (Nr +hk ' n{x \^)dA dx = 0, t > 0, 99 G S(R), if > 0. (3.19) 



T 



Summarizing, to obtain Theorems 2.2 and 2.4 we prove convergences of (|3.14p - (|3,16p and (|3.19p . 
In each case these proofs require some non-trivial work. 

3.3 Proof of Theorem 2.2(a) 

Following the scheme we show that 



lim I(T;k,n) = pi(0) / / |r — r'\ 1 l a dr'dr \ ip(x)dx / ip(x)da 

JO JO JR JR 



T— > 

lim H(T;k,n,m) = lim M(T;k,n;£,m) 

T->oo T-^oo 



(3.20) 



= Pl(0) / / ( r + r'y 1/a dr'dr [ (p(x)dx [ ip{x)dx. (3.21) 
Jo Jo Jr Jr 

(see (I344"]1 - (I3~T61) N 1. It is easy to see that by (|3TT3|) . and ([13]), this yields (fHTTUD . 

Let 77 denote the standard a-stable Levy process in IR, As we consider the model without branching, 
we have, for r > r', 

= Eif(x + h k ,n(x) + 7] r )lp(x + h ktn (x) + T] r >) 

= T r '(lpT r -r"f)(x + h ktn (x)). (3.22) 

Putting this into (|3.14|) and omitting the subscripts k, n, we obtain 

I{T) = h(T) + I 2 (T), (3.23) 
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where 

h(T) = —[ I li r>r >\[ p r >{x + h(x) -y)ip(y)p r _ r ,(y - z)tp(z)dzdydxdr'dr, (3.24) 
JO JO Jm 3 

h(T) = -=2 / / l{r<r'}/ p r (x + h(x) - y)(p(y)p r ,^ r {y - z)4>(z)dzdydxdr'dr. (3.25) 
Jo Jo Jr 3 

In Ii we substitute r = r/T, r' = r'/T, use ()2.5p and (13. ip . and then we substitute x = T~ l l a (x—y), 
arriving at 

h(T) = f f\ {r>r , } ! p r ,{x + T- l l a h{T l l a x + y)W(y) 
Jo Jo Jr 3 

■p r ^ r i(T~ 1 ' a (y — z))ip(z)dzdydxdr' 'dr. 
By (|3.17p and (|3.ip . the expression under the integrals converges pointwise, as T — > oo, to 

Pr'{x)p r -r'(.0)^{y)ip(z) = p r ,(x)\r - r'^^px (Q)ij){y)ip (z), 

and by (|3.3p . for T > 1, it is bounded by g r >{x)\r — r'\~ 1 ^ a pi(0)ip(y)ip(z), which is integrable over 
[0, t] x [0, s] x M 3 , since a > 1. .Z2CO can be treated analogously, hence by (|3.23p we obtain (|3.20|) . 
Next we take E. In f|3. 15|) we substitute r = r/T,r' = r'/T, and by (|2.5p we have 



E(T; k, n, m) 



t r s 



JO 



T 1/a E / / / p Tr (x + hk, n (x) -y)prr'(x + h k>m (x) - z)if(y)tj;(z)dxdydzdr'dr. 



We use (13. ip and substitute x = T 1 l OL {x — z), obtaining 

I(T;k,n,m) = E I f / p r (x + T~V a {z - y + h Kn {T l l a x + z)) 
Jo Jo Jr 3 

■p r >{x + T- 1/a h k , m (T 1/a x + z))tp{y)4){z)dxdydzdr' dr. 

The integrand converges pointwise to p r (x)p r / (x)(p(y)ijj(z) , and ()3.17|) implies that for T > 1 , it is 
bounded by r~ 1 / a pi(0)g r '(x)if(y)ip(z) (see (|3.3p ). which is integrable. As J R p r (x)p r /(x)dx = (r + 
r')~ 1/cl pi(0), we obtain the limit (pHT]) for E. 

Note that in this argument the only property of h k , n we have used is (I3.17p . therefore it is immedi- 
ately seen that the limit of M can be obtained in the same way (see (|3.16p ). This completes the proof 
of (l3~TUjl . 

As explained in the scheme, to finish the proof it suffices to show (|3.19p . 
The expression under limy— >oo 

sup n k in (|3.19p . similarly as in (I3.22j) . can be written as 



Tt r Tt rTt 

T r ((pT r '-r( l P'Tr"-r' i p))(x + hk, n (x))dr" ' dr' ' drdx 



3\_ 

Ft j Jr Jr 1 



6 r r Tt r Tt r Tt f 
— "itf / / / / 7 r ( l P%'(f%"f))(x + hk, n (x))dr dr drdx 
F T JrJo Jo Jo 

< CJ(T), (3.26) 

where 



1 r rTt rTt rTt 

J(T) = -3 / / / / Tr(foTr>(fcT r »fo))(x)dr"dr'drdx, (3.27) 

?T Jo Jo Jo 
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and 4>2 is given by (I3.4p . The last inequality in (13.261) is of the same type as (j3.6|) . and can be obtained 
by an analogous argument using f)3.5j) and (|3.17p . Hence, for f)3. 19|) it is enough to show that 

lim J(T) = 0. (3.28) 

T— s-oo 

Note that in this argument we have not used the assumption on a. 

After obvious substitutions, using (|2.5p and the invariance of Lebesgue measure for 7t, we have 

J(T) = T 3 / 2a f f f [ Mx)PTr>(x ~ y)My)PTr»(y - z)fo(z)dzdydxdr"dr'dr. 
Jo Jo Jo Jr 3 

By the self-similarity property (13.ip . 

< T 3/2a-2/* p 3^ t n\-V<x ds \ (J 

since 1 < a. □ 
3.4 Proof of Theorem 2.2 (b),(c) 

According to the scheme, we prove A3. lOf) . Fix s < t. (I3.14p can be written as 

I(T;k,n) = h(T;k,n) + I 2 (T;k,n) + I 3 (T;k,n), (3.29) 



3 

ix ) — > as T — > oo, 



where 



j rTs rTs r 

Ii(T;k,n) = -p^E III %■{,¥>%•' -r^)(x + hk ifl {x))dxdr' 'dr, (3.30) 
F T Jo J r Jr 

j rTs rr r 

l2(T;k,n) = — I II T r '(i^Tr-r ,( p)(x + hk >n (x)) dxdr'dr, (3.31) 
Ft Jo Jo Jr 



T JO JO 
i r-Tt r Ts 

t JTs Jo 



I 3 (T;k,n) = —^E \ \ \ ...dxdr'dr, (3.32) 



where ... is the same integrand as in (|3.3ip . 

To compute the limit of I\ we fix an arbitrary < e < 1, substitute r = r/T, and we write 

h(T;k,n) =l' 1 (T;k,n)+l'{(T;k,n), (3.33) 

rp pes rTs p 

I[(T;k,n) = —^E III Trr^pTr' -Trtp){x + h ktn (x))dxdr'dr, (3.34) 
£ T Jo Jtt Jr 



l"(T;k,n) = —^E / / / ...dxdr'dr. (3.35) 
F T Jes Jtt Jr 

Consider case (c): a < 1. Applying a version of (|3.6p we have 

rp res rTs r 

l[(T;k,n) < —^C / / / T T r{<p2T r >(t)2){x)dxdr'dr 
£ t Jo Jo Jr 

= -t^Ces ( j <p2(x)T r '(f>2(x)dxdr' (3.36) 
F T Jo Jr 

< Cis [ fa{x)Gfa{x)dx = C 2 e, (3.37) 
Jr 
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by H32J and §ZB. Next, 

rp r-s rT(s-r) 

7 T Jes JO 



l"(T;k,n) = I I I „PTr(x + h^ n (x) - y)tp(y)%>ip(y)dydxdr'dr (3.38) 



fT(s-r) 

£ / / Pr(.x + T-V«h k , n (T 1 ' a x + y))tp(y) I T r ^(y)dr'dydxdr, 



where we have applied (|3.ip and the substitution x = T- l / a {x-y). 
It is now easy to see that 

lim I"(T; k, n) = (s- es) [ ip(y)G^(y)dy, (3.39) 
T -^°° Jr 

since the passage of the limit under the integrals is justified by f|3.3[) . Note that the function g r {x) is 
not integrable in r in a neighborhood of for a < 1; that is why we had to consider the interval [es, s]. 
Prom ([3T33]) . (|33Tj) and flES}, we infer that 



lim h(T;k,n) = s / <p(y)G*l>(y)dy. (3.40) 

Note that ()3.31|) is essentially the same as ()3.30p . only the roles of <p and ip are interchanged. Therefore, 
by symmetry of G, we have 

lim I 2 (T;k,n) = s [ ip(y)Gi{j(y)dy. (3.41) 

Passing to /3 defined by (3.32), we first estimate it similarly as I\ (see (|3.36p ). then we change the 
order of integration drdr', and substitute r' = r'/T and f = r — Tr', obtaining 



ps poo p 

Is(T;k,n)<C / / / 4>2{x)T r 4>2{x)dxdrdr' . 

Jo Jr(s-r') Jr 

Now, it is clear that 

lim Iz(T\ k, n) = 0. 

T->oo 

hence, by (j339|) . p^Ojl and (g3H) , we have 

lim I(T; k, n) = 2(s A t) [ ip{x)Gil>{x)dx. (3.42) 
Next, by ([BTToD . (^TdT) . using (pPTjl and ([22]), we obtain 



C # /* s f 

E(T;k,n,m) + M(T;k,n;l,m)< — ^ / / / 7r4 > 2{x)l~ r "l ) 2{x)dxdr' dr 

Ft Jo Jo Jr 

4>2(x)% +r i<f>2{x)dxdr'dr (3.43) 



T 



"t j o Jo 

pt poo p 

< C / Mx)T r >Mx)dxdr'dr -> 0, (3.44) 

Jo JTr Jr 

since, by (|3.7p . G(/>2 is bounded. This together with (|3.42p and ([3.13P proves (|3.10p for a < 1. 
Now consider case (b): a = 1. We also have (I3.36p . hence we get 

l[{T;k,n) < Ce, 
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because for such a, 



i r 

sup sup- — — / %4>2 (x) dr < oo, (3.45) 

T>2x£R K>gi J 

see (3.46) in [5]. It is easy to see that the limit of Ii(T;k,n) given by (|3.38p is the same as the limit 
of 



1 



logT 



s-l/T rT(s-r) 



es J 1 



El I I pTr(x + h k , n (x) - y)(fi(y)p r '(y - z)ip(z)dzdydxdr'dr. 



Using (|3.ip twice and substituting r' = log r 1 / log T and x = T 1 {x — y), this expression is equal to 

/•s-l/T Aog T(s— r)/ log T /• 

E / p r (x + T~ 1 h(Tx + y))(p(y)pi(T~~ r (y — z))ifj(z)dzdydxdr'dr 



es JO 



— > s(l-e)pi(0) / / ip(z)dz, 

by (|3.3|) . Thus, we have shown that 



lim I\ 

T->oo 



(T;k,n) = (sAi)pi(O) / ^(x)cte / ^(y)dy = lim I 2 (T; fc, n), (3.46) 



where the last equality follows by symmetry (see (|3.3ip ). 

To estimate ^3 given by (3.32), we again use an inequality of the type (|3.6p . obtaining 

I 3 (T;k,n) < — — — / / / p r >(x - y)4> 2 (y)p r - r '(y - z)4> 2 {z)dzdydxdr' dr 



TlogT J Ts 







C ( f \ 2 f Tt [ Ts 

< — — ^Pi(O) / 4> 2 {y)dy) I I (r - r'Y x dr'dr -> as T 00. 
TlogT \J R ) J Ts Jo 



Similarly, by (|3.43p we have 



E(T;k,n,m) + M(T;k,n;£,m) 



< -^—pi(0) f [ {r + r')- x dr'dr( f 4> 2 (x)dx 



logT 



Jo 



— > as T — > 00. 

This, together with ([3^61) . (I3T29D and (I3T3D prove (I3T0D . 

It remains to show (|3.19|) . To this end, by (|3.26p and (|3.27|) . it suffices to observe that (|3.28|) holds. 
Indeed, for a < 1 we use boundedness of G(p 2 , and for a = 1 we employ (|3.45p . The proof of the 
theorem is complete. □ 

3.5 Proof of Theorem 2.4 

We recall first the following formula for the second moments of critical binary branching systems with 
branching rate V: 

E{N?,<p){N?,1>) 

rr 

= T r (vT r i-rlf>) + V / T u {{T r -uV){Tr>-ulp))du, T < r' , X G R, (3.47) 







which is obtained using, e.g. Lemma 3.1 in [25], and the Markov property. 
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Fix s <t. By (13.47H . expression (I3.14|) can be written as 

J(T; k, n) = Vh(T; k, n) + I 2 (T; k, n), (3.48) 



where 



1 f s f f s f 

h(T;k,n) = -^E / / Tu{{T r -uV){T r ,- u ^)){x + h k>n {x))dxdr> 'drdu, (3.49) 

JO Ju J^- 

I 2 (T;k,n) = -^Ej f f {T r {^pT r ' - r ^){x + h kjn (x))l {r < r n 

?T Jo Jo Jr 



+T r i(i/j%- r 'ip)(x + hk : n( x ))^{ r ' <r})dxdr' dr. 



Comparing 1 2 to (|3.29p - (|3.32p . we see that it has the same form as I(T;k,n) in the non-branching 
case with a < 1, only the norming Ft may be different. It is now immediately seen that by (13.42|) we 
have 

T lim 7 2 (r;fc,n) = | 2 ( 8 At) f R X)G^(x)dx, if l<l<\. (3 ' 50) 
Similarly, using the result for the non-branching system (see f|3.44j) ) . we obtain 

lim E(T;k,n,m) = lim M(T;k,n;£,m) = 0. (3.51) 

T— >oo T->oo 

Hence, to prove (|3.10p it suffices to calculate the limit of I\(T; k,n) given by (|3.49p . 
Consider the case 1/2 < a < 1. Making the substitutions f = (r — u)/T,r' = (r' — u)/T,u = u/T, 
using (|3.ip and (|2.12p . and then putting x' = T~^' a x,y = T~ l < a y, we have 

Ii(T; k,n) 

ps pt—u rs—u p 

= T^ a E / / / Pu(T- 1/a (x + h k , n (x)-y))p r (T-V a (y-z)) 
Jo Jo Jo Jr 4 

■p r i(T~ l l a (y — w))ip(z)ip(w)dwdzdydxdr' drdu (3.52) 

f-s rt—u i-s—u 



ps pt—u ps—u 

E / R T ,k,n(u,r,r')dr'drdu, 

Jo Jo Jo 



where 



R T ,k,n(u,r,r') = [ p u (x-y + T- l l a h k , n {T l l a x))p T {y-T- l l a z) (3.53) 
Jr 4 

'Pr'(y — T~ 1 l a w)ip(z)'4)(w)dwdzdydx. 
Using (|3.2p . (p < Ccfim, (|3.5p and (|3.ip . it is easy to see that 

lim R Tk n (u,r,r') = / p u (x - y)p r (y)p r >(y)(p(z)i/j(w)dwdzdydx 

T^oo J R 4 



= (r + r') _i/ >i(0) / <p(z)dz / ip(w)dw. 

Jr Jr 

By (I332|) . (I33TT) . (IQ81) . (l330|) . fl343"l) and ([23]), it is clear that in order to obtain (ETTOj) it remains 
to justify the passage to the limit under the integral in ()3.52p . To this end, in (|3.53p we substitute 
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y = y-T- l / a h Kn {T l l a x),~z = z- h k , n (TV a x), w = w - h^T^x), use ^(z + ^^T 1 /^)) < Cfa(z) 
(see (|3.17p and (|3.5p ). and the same for -ip, obtaining 



RT,k,n( U ' r i r ') ^ C Pr+r'{T 1/a {z - W))(j> 2 {z)(t)2{w)dwdz 

Jk 2 

< C(r + r')- 1/a pi(0) H ' <fo{z)d^j . 

The last expression is integrable on the set of integration in (|3.52p . since a > 1/2. 

Now consider the case < a < 1/2. After obvious substitutions, and using = T,I\ given by 
(|3.49p can be written as 

ps rT(r—u) pT(s-u) p 

Ii(T;k,n) = E / / Ttu{{T t <p){T t 'iP)){x + hk, n (x))dxdr' drdu 



by (pTfjl . (1331) and (ETTTD . hence 



since a < 1/2. 



Similarly as in (I3.38P we obtain 

lim l"(T;k,n) = (s - es) / G(p(y)Gijj(y)dy, 

hence, by ([337jl and (j3"lg]) . 



lim 

T->oo 



h(T;k,n) = (sAt) f G^(y)G^(y)dy. 

JR 



JO JO JO JR 

= I{(T;k,n) + l'{(T;k,n), (3.54) 
where, for any < e < 1, 

/■es 

l[(T;k,n) = E (3.55) 

JO 

l"(T;k,n) = E f (3.56) 

(cf. (13331 - (jS35]» . We have 

/{(T; k,n) < E / p Tu (y)Gip(x + /i fc .„(x) - y)Gip(x + h k , n (x) - y)dydxdu 

Jo Jm 2 

< C / PTu(y)—— ; ion n ^ d vdxdu, 

Jo J m 2 1+ s -y 2U-") 



I[(T;k,n) < Cte, (3.57) 



This, together with (l348|) . (l330|) . (j33Tjl . (j3X3jl and (^T3jl . proves (I3T0]) . 

It remains to show ([3.10P in the case a = 1/2, which amounts to proving that 

lim Ii(T; k,n) = (s A t)— [ ip(x)dx [ ip(x)dx. 

We decompose I\ as in (|3.54p - (|3.56p . The integral (T/F^)E J Q £S ... is easy to estimate. To obtain 
the limit of (T/F^)E J^ g ... we combine the argument on pp. 1354-1356 of [10] and the method for 
getting rid of hk, n (x), which we have used several times above. We skip the cumbersome details. 



16 



We have proved (I3.10P in all the cases. According to the scheme, to finish the proof it remains to 
show (|3.19p . To this and we define the function 

v e (x,t) = I- Eexpl-e J (N?,<p)dr\ , 9 >0,x £R,t>0. 

It is known that the Feynman-Kac formula implies that vq satisfies the non-linear equation 

ve(x,t) = J Tt-u \ 0<p{-)(\ - v e (-,u)) - ^vq(-,u) \ {x)du 

(see e.g. [18], or the space-time approach used in [U[5]). Hence, by a similar argument as in (3.45)- 
(3.47) of [5], we obtain 

E (£\N^ +hk ' n(x \ V )dA = E-^v e (x + h kin (x),Tt)\e= 
= 6E J Trt-r yf J T- u yp J T s <pds\ di?j (x + h k<n {x))dr 
+ WE J Trt-r {^P J %-u(^J T s ipdsj dv^j (x + h kiU (x))dr 
+ 6VE J Trt-r (^J Tpdv J %-u (p J T s <pds \ dv^j (x + h k ^ n (x))dr 
+ 3V 2 E J Trt-r yj T v (pdv J T r - U {^J Tpds^j duj (x + h k>n (x))dr. 
Without loss of generality we may assume t = 1. 

Repeating the argument as in (|3.5p . (|3.6p . and using the invariance of Lebesgue measure for Tt, it 
is not hard to see that in order to prove (I3.19P it suffices to show that 



where 



lim Ji(T) = 0, i = 1,2, 3, 4, (3.58) 

T— >oo 



J\(T) = -^3 J j (j) 2 {x) j Tr-u ($2 j TsfadsJ (x)dudrdx, (3.59) 
■h{T) = J J <p2(x) J %-u {^J Tfods^j (x)dudrdx, (3.60) 



Jz{T) = — 3 / / / T v 4>2{x)dv l T-u ( 4*2 / Tsfads ) (x)dudrdx, (3.61) 

■'K J J J \ JO 

T rr rr / ru \ 2 



1 

7 T JM JO ■/() •/() 

Assume 1/2 < a < 1. Denote 



J 4 (T) = -^3 / / / TMx)dv j T-u ( / Tfods j (x)dudrdx. (3.62) 



/(*) = f 1 pr(x)dr, 4> T (x) = T 1 / a MT- 1/a x). 
Jo 

We have 

||HK = ||^2||i, ||0r||2 = T 1 / 2a ||</> 2 || 2 , (3.63) 
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(|| • Hp denotes the norm in L P (M)) and 

||/||2<oo, (3.64) 

since a > 1/2. The proofs of (|3.58p use (|3.63p and (|3.64|) together with the Schwarz and Young 
inequalities. For brevity, we show (13.58P for J% and J4 only. 
By (j339|) and (|215]l . 



MT) < T- 7 / 2+3 / 2a [ <f> 2 (x) / T u (<h [ Tsfcds) (x)dudx 

JR JO \ JO J 

= T- 3 / 2 - 1 / 2 * I Mx)f(x-y)My)f{y-z)4> T (z)dzdydx, 



by obvious substitutions and (|3.ip . Hence 



j l{ T) < T- 3 / 2 -V 2a \\fo\\ 2 \\f*(Mf*t T ))\\ 2 

< r-^-v^n^iiiii^ii,! 

< CT- 3 ' 2+l l 2a -> 0, 



by (I3T63D . dSMD , and since a > 1/2. 

Similarly, by (|3.62p and (I2.12P we have 



J 4 (T) < T - 7 / 2+3 / 2a [ [ T %Mx)dv [ T T u ( [ T Tsfods) {x)dudx 

Jr Jo Jo \Jo J 

= TV 2 -V 2 «\\U*]>T)U*U*j>T?)\\l 

< r 1 / 2 - 1 /2«|| / ||3||^ 2 ||3^ 0) 

since a < 1. The remaining limits in (13.58P are obtained similarly. 
Now assume a < 1/2. By (I3.7|) we have 

G(f>2 G L°°(K.) n L 2 (R) and G((G^ 2 ) 2 ) e L°°(R) n L 2 (M). (3.65) 

These properties easily imply (I3.58p . For example, 

Ji(T) < -jL / <h(x)G(<hG<h){x)dx < 4=, (3.66) 

and 

J 4 (T) < I3 / G0 2 (x)G((G0 2 ) 2 )(x)dx < r ' 



Finally, assume a = 1/2. In this case we use G(f>2 £ L°°(M) and 

_ / 1 / ' 

log: 



-*—:[([ T u tp(x)du) dx<C for T > 2, (3.67) 
%T J R \J ) 



see (3.33) in [7J. 

It is clear that (|3.66p also holds (recall that Ft = \/T logT) and (|3.58j) for J 2 and J3 follows easily 
from (|51T7jt . 

We turn to J4, which requires more work. By (I3.62p . the Schwarz inequality and (I3.67P we have 

J 4 (T)< ] (3-68) 
VT log T 
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where 



R{T) = J(JTu (J Tsfods^j (x)du\ dx 



j 2 J J Pu+u'{y ~ z)du'duij T s (j>2{y)ds\ (J T s > folds' J dzdy 
R 1 {T) + R 2 (T), 



where 



Rl (T) = Iff--. 

JR 2 JO JO 

R2(T) = [ f f ... 
Jr 2 Ji Jo 



We have 



'/' /■/' / /- / /■:/' \ 2 \ 2 

T s <j)2(y)ds I dy 



R2(T) < pi(0) f [ (u + u')- 2 dudu' \ [ ( [ 

J l JO \JR \J0 

< C(logT) 3 , 

by dSHTJ). Using 



Pu+u'{y - z)dudu < 



Jo 



\y-z\V* 



we obtain 

R X {T) < R[(T) + R'{{T), 

where 

1 



R[(T) = Cj J J (G<j> 2 (y)f(Gct> 2 (z)) 2 dydz, 

\y-z\<l 

R'((T) = C f f (j\sMy)ds^J n\ s ><f>2(z)ds^ dydz. 

\y-z\>l 

For R'i(T) we use (]3.7p and (I3.5p . obtaining 

R'l(T) <Cx f 1 dy [ -^dx < oo, 

Jr (i + \y\) J\x\<i m 1 ' 2 

and by (^671) . 

R'{{T)<C 2 (\ogT) 2 . 
Putting these estimates intro (|3.68p we arrive at 

which completes the proof of the theorem. □ 
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3.6 Proof of Proposition 2.6 



Since the convergence of finite-dimensional distributions has been proved already, by virtue of the 
Mitoma theorem [30] it remains to show tightness of (Xt,<P),T > 1, in C([0, t],M) for any fixed 
cp £ S(M), if > 0. To this end we prove 

E((X T (t), ip) - (X T (s), ip}) 2 < C(t - s) a , s<t<r, (3.69) 

for some a > 1. 

By ([3SD and (|3J| we have 

E((X T (t),<p) - {X T (s)^)f 

rj+l fTt fTt 



i p+i fit fit 
^ El?*/. / / E(N?,<P)(N^<p)dr'drdx 

jp.7, T J 3 JTs JTs 



1 



Tt r Tt 

E(N^,(p)(N^,(p)dr'drdx. 



T JTs JTs 

(P) (P) 9 (P) 

The latter expression is equal to E(Xj-, (t) — Xj, (s),ip) for Xj, corresponding to the system 
starting from the standard Poisson measure (with or without branching), and we know that (I3.69P is 
satisfied for X^ (see subsection 3.1 in [1]). 

Remark Recall that Proposition 2.6 refers to a special simple choice of v. For a general v £ A4, the 
proof of ()3.69p is similar but slightly more involved. One has to estimate an extra term and use an 
inequality of the type (|3.6[) . 
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